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ABSTRACT 

We explore the growth of <10'' Mq black holes that reside at the cen- 
ters of spiral and field dwarf galaxies in a Local Group type of environment. 
We use merger trees from a cosmological N-body simulation known as Via 
Lactea II (VL-2) as a framework to test two merger-driven semi-analytic 
recipes for black hole growth that include dynamical friction, tidal strip- 
ping, and gravitational wave recoil in over 20,000 merger tree realizations. 
First, we apply a Fundamental Plane limited (FPL) model to the growth 
of Sgr A*, which drives the central black hole to a maximum mass lim- 
ited by the Black Hole Fundamental Plane after every merger. Next, we 
present a new model that allows for low-level Prolonged Gas Accretion 
(PGA) during the merger. We find that both models can generate a Sgr 
A* mass black hole. We predict a population of massive black holes in local 
field dwarf galaxies - if the VL-2 simulation is representative of the growth 
of the Local Group, we predict up to 35 massive black holes (< 10^ Mq) 
in Local Group field dwarfs. We also predict that hundreds of < 10^ Mq 
black holes fail to merge, and instead populate the Milky Way halo, with 
the most massive of them at roughly the virial radius. In addition, we find 
that there may be hundreds of massive black holes ejected from their hosts 
into the nearby intergalactic medium due to gravitational wave recoil. We 
discuss how the black hole population in the Local Group field dwarfs may 
help to constrain the growth mechanism for Sgr A*. 
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Local Group, AGN feedback, supermassive black holes, dark matter halos, 
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1 INTRODUCTION 

Supermassive black holes, with masses of IO^Mq < 
M < 10^'' Mq, are widely believed to dwell at the cen- 
ters of elliptical galaxies and spiral bulges (e.g. Kor- 
mendy & Richstone 1995); the best known example is 
observed at the center of the Milky Way, with a mass 
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MsMBH = 4.2x10^^ M0 (Ghez et al. 2008). There is 
abundant evidence that when a SMBH is in place, it 
transforms the structure and evolution of the galaxy, 
from powering active galactic nuclei at high redshifts 
(Greenstein & Matthews 1963, Rees 1984, Alexander 
et al. 2005, Fan 2005), to regulating star formation 
throughout the galaxy (Di Matteo et al. 2005, Croton 
et al. 2006, Cox et al. 2008), to scouring the galactic 
nucleus of stars during SMBH mergers (Ebisuzaki et 
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al. 1991, Quinlan 1996, Makino 1997, Milosavljevic & 
Merritt 2001, Volonteri et al. 2003). 

This deep connection between the evolution of 
SMBHs and galaxies is perhaps best encapsulated in a 
remarkable correlation between the SMBH mass and 
the velocity dispersion of the host spheroid (Gebhardt 
et al. 2000, Ferrarese & Merritt 2000, Tremaine et al. 

2002, Marconi & Hunt 2003, Hopkins et al. 2007). The 
dispersion in the black hole fundamental plane (Hop- 
kins et al. 2007) points to an intrinsically tight corre- 
lation, at least for a sample of nearby bright spiral and 
elliptical galaxies with clear dynamical SMBH signa- 
tures. However, on a smaller mass scale, in systems 
that have a mass comparable to the Milky Way mass 
or smaller, central SMBHs may become less common 
as bulges become less common (Ferrarese et al. 2006), 
and in some instances they disappear entirely, such as 
in the case of M33 (Merritt et al. 2001, Gebhardt ct 
al. 2001), and NGC 205 (Valluri et al. 2005). Many 
of these bulgeless stellar systems host a nuclear star 
cluster instead. Nuclear star clusters are found in late- 
type spirals (Bokcr ct al. 2002), and dwarf elliptical 
galaxies (Cote et al. 2006). This changeover may be 
the result of a competition between the SMBH and 
nuclear star cluster for the same gas reservoir. Nayak- 
shin et al.(2009) shows that in massive galaxies with 
a spheroid velocity dispersion of cr > 150 kms^^, gas 
accretion onto the black hole dominates and a SMBH 
forms, while in galaxies with lower velocity disper- 
sions, star formation generates a nuclear cluster. Note 
that such a star cluster does not exclude presence of 
an underweight SMBH (Nayakshin et al. 2009). 

Several observational and theoretical studies have 
linked the SMBH mass to the mass of the host dark 
matter halo (Ferrarese 2002, Baes et al. 2003, Shankar 
et al. 2006). This relation is a boon to theorists be- 
cause many of the leading explanations of SMBH birth 
and growth are driven by hierarchical structure for- 
mation (Hopkins ct al. 2005, Wyithe & Locb 2005, 
Granato et al. 2001, Menou et al. 2001, Adams et al. 
2001, Monaco et al. 2000, Silk & Rees 1998, Haehnelt 
et al. 1998, Haehnelt & KaufFmann 2000, Cattaneo et 
al. 1999, Loeb & Rasio 1994), and are therefore tied 
to the mass of the dark matter halo. 

In the current picture of SMBH assembly, the 
black hole begins life as a low mass "seed" black hole 
at high rcdshift. It's not clear, though, when exactly 
these BH seeds emerge or what mass they have at 
birth. SMBH seeds may have been spawned from the 
accretion of low angular momentum gas in a dark mat- 
ter halo (Koushiappas et al. 2004, Bromm & Loeb 

2003, 2004), the coalescence of many seed black holes 
within a halo (Begelman & Rees 1978, Islam et al. 
2004), or from an IMBH formed, perhaps, by runaway 
stellar collisions (Portegies Zwart et al. 2004, Miller 
& Colbert 2004, van der Marel 2004) or they could 
even be primordial (Mack et al. 2007). However, the 
most likely candidates for SMBH seeds are the rem- 
nants that form from the first generation of stars sit- 



ting deep within dark matter halos (Madau & Rees 
2001, Hcgcr et al. 2003, Volonteri ct al. 2003, Islam 
et al. 2003, Wise & Abel 2005) - so called Popula- 
tion 111 stars. With masses less than roughly lO"' Mq, 
these relic seeds are predicted to lie near the centers 
of dark matter halos at high redshifts (Bromm et al. 

1999, Abel et al. 2000, 2002). Structure formation dic- 
tates that dark matter halos form in the early universe 
and hierarchically merge into larger bound objects, 
so naturally as dark matter halos merge, seed black 
holes sink to the center through dynamical friction 
and eventually coalesce. 

Gas accretion is thought to play a critical role in 
fueling the early stages of black hole growth (David et 
al. 1987, Kauffmann & Haehnelt 2000, Merloni 2004), 
and this may explain the tightness of the Mbh — cr 
relation (Burkert & Silk 2001, Haehnelt & Kauffmann 

2000, Di Matteo et al. 2005, Kazantzidis et al. 2005, 
Robertson et al. 2006). Since high redshift galajcies are 
thought to be especially gas-rich, each merger brings a 
fresh supply of gas to the center of the galaxy, and new 
fuel to the growing supermassive black hole (Mihos 
& Hernquist 1994, Di Matteo et al. 2003). From a 
combination of gas accretion and binary black hole 
coalescence, it is thought that these Pop Ill-generated 
seeds may form the SMBHs we observe today (Soltan 
1982, Schneider et al. 2002). 

During a galaxy merger, each black hole sinks to 

the center of the new galaxy potential due to dynami- 
cal friction and eventually becomes bound as a binary 
(Kazantzidis et al. 2005; Escala ct al. 2005). Dynam- 
ical friction then continues to shrink the orbit until 
the binary is hard (i.e, the separation between each 
black hole, aBSH , is such that the system tends to lose 
energy during stellar encounters) (Heggie et al. 2007). 
Thereafter, further decay is mediated by 3-body scat- 
tering with the ambient stellar background until the 
binary becomes so close that the orbit can lose energy 
via gravitational radiation. In studies of static, spheri- 
cal potentials, it may be difficult for stellar encounters 
alone to cause the binary to transition between the 3- 
body scattering phase and the gravitational radiation 
regime (Milosavljevic & Merritt 2003). However, in 
gas-rich or non-spherical systems, the binary rapidly 
hardens and coalesces into one black hole, emitting 
copious gravitational radiation in the process (Mayer 
et al. 2007, Kazantzidis et al. 2005, Berczik et al. 
2006, Sigurdsson 2003, Holley-Bockelmann & Sigurds- 
son 2006). 

Unless the black hole binaries stall at the final 
parsec, the longest timescale governing the coalescence 
of two black holes occurs when the host galaxies them- 
selves are still merging. Here, the dynamical evolution 
of two merging galaxies is driven by the combined ef- 
fect of dynamical friction, which brings the less mas- 
sive galaxy, or satellite, to the center of the larger 
halo, or primary, and tidal interaction which strips 
mass from the satellite and further delays the merger 
(e.g. Richstone 1976, Aguilar & White 1986, Holley- 
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Bockclmann & Richstonc 1999, Taffoni et al. 2003). If 
the dynamical friction timescale is longer than a Hub- 
ble time, the black holes carried by their host galaxies 
will not sink close enough to form a binary. 

Sijacki et al. 2007 and Di Matteo et al. 2008 
have performed state-of-the-art high resolution hydro- 
dynamic simulations of cosmological structure forma- 
tion, following the growth of high mass SMBHs at the 
centers of massive elliptical gala^xies and clusters of 
galaxies. Their research was followed by similar semi- 
analytic work that incorporates a full treatment of 
dark matter dynamics, radiative gas cooling, star for- 
mation and energy feedback processes (Somerville et 
al. 2008). In this very elegant approach, the SMBH 
accretes gas through a quasar mode - nearly Edding- 
ton rate accretion following a galaxy merger (Croton 
et al. 2006) - and a radio mode - Bondi-Hoyle ac- 
cretion associated with relativistic jets (Somerville et 
al. 2008). Both modes produce feedback that heats the 
surrounding gas. In this model, the feedback stops the 
accretion and locks the growth of the SMBHs to the 
fundamental plane. At the same time, the feedback 
also quenches the star formation, which explains the 
observed shallow metallicity, stellar density and en- 
tropy profiles. However, it has recently been suggested 
that the implementation and importance of AGN feed- 
back may need to be reexamined (Ostriker et al. 2010). 
Also, radiation fields and winds produced by massive 
stars may provide the dominant feedback (Hopkins et 
al. 2010). In fact, it is worth mentioning that AGN 
feedback is just one of many feedback processes oc- 
curring in galaxy centers; supernovae, star formation, 
and galaxy mergers produce feedback as well (Sinha & 
HoUey-Bockelmann 2010). It is not clear which feed- 
back process contributes the most to the evolution of 
a given galaxy. 

As well-developed as the current effort to under- 
stand SMBH growth is, most work has focused on 
growing the most massive SMBHs, found predomi- 
nantly in dense environments. In a Local Group en- 
vironment though, the galaxy morphology, dynamics, 
and star formation history arc all dramatically differ- 
ent. Smaller stellar systems such as disk and dwarf 
galaxies, have SMBHs with lower masses (< 10'' M©) 
or no SMBHs at all for bulgeless galaxies like M33. In 
cosmological simulations, AGN feedback models cre- 
ate significant uncertainties for black hole growth in 
these lower mass halos (Booth & Schaye 2009). In 
addition, SMBHs in the Local Group are currently 
subluminous, implying a different relationship with 
gas consumption, at least in the present epoch, than 
the more massive end of the SMBH mass spectrum. 
Nonetheless, mergers are still expected to provide a 
key mechanism to channel gas to the galactic center. 
While "cold mode" gas accretion (Keres et al. 2005, 
Dekel et al. 2009) can channel large amounts of cold 
gas (< 10* K ) to the galactic disk, this gas may not 
fuel the SMBH. Indeed, high resolution N-body-SPH 
simulations of Milky Way formation show that most 



of the gas that is accreted originates from major merg- 
ers; even during the "radio mode" , where the gas is ac- 
creted by a Bondi - Hoyle mechanism, the gas reservoir 
is stocked not from the cold mode, but from the major 
mergers the Milky Way finished billions of years ago 
(J. Bellovary, private communication). Hence, galaxy 
mergers appear to drive SMBH growth in a wide range 
of masses. 

One problem that can preferentially plague low 
mass galaxies is gravitational wave recoil. Recent cal- 
culations of binary black hole mergers indicate that 
gravitational wave recoil can kick a newly merged 
black hole with a speed as large as ~4000km s~^ (Her- 
rmann et al. 2007, Gonzalez et al. 2007a, 2007b, Kop- 
pitz et al. 2007, Campanelli et al. 2007, Schnittman & 
Buonanno 2007). The magnitude of the gravitational 
wave recoil imparted during any given merger depends 
on the mass ratio of each black hole, the spin magni- 
tude and orientation with respect to the binary orbital 
plane, as well as the eccentricity of the orbit (Campan- 
elli ot al. 2007, Schnittman & Buonanno 2007, Baker 
et al. 2007). While certain spin alignment mechanisms 
(Bogdanovic et al. 2007, Dotti et al. 2010) could allow 
massive galaxies to retain their SMBHs, even mod- 
erate kicks can eject a growing SMBH from dwarf 
or high redshift galaxies (Merritt et al. 2004, Micic 
et al. 2006, Volonteri 2007, Schnittman 2007, Sesana 
2007, Volonteri et al. 2010). There is tentative observa- 
tions that suggest that gravitational wave recoil can 
indeed eject a SMBH from its host (Komossa et al. 
2008). Given the greater vulnerability to recoil, the 
questionable importance of AGN feedback, and lack 
of SMBH in some galaxies, one can argue that SMBH 
growth may be substantially different for these light- 
est SMBHs in Local Group type of environments. 

In this paper we use Via Lactea 2 (VL-2) evo- 
lutionary tracks to begin with an N-body merger tree 
that mimics the assembly of the Milky Way halo (Die- 
mand et al. 2007,2008). Upon this numerical merger 
tree we paint semi-analytic recipes for the galaxy and 
black hole growth. We follow two different merger 
driven recipes. The first approach is inspired by re- 
cent semi-analytic and rmmerical work on black hole 
growth (Somerville et al. 2008, Croton et al. 2006, 
Hopkins et al. 2005, 2007). Here, SMBHs grow slowly 
through low-level Bondi-Hoylo gas accretion tempered 
by AGN feedback - the radio mode; this growth is 
punctuated by rapid Eddington gas accretion trig- 
gered by a major merger - the quasar mode. This 
approach has been used to explain SMBH growth in 
massive stellar systems but not the growth of massive 
black holes in the Local Group. Here, SMBH growth 
is limited by the black hole fundamental plane. Hence, 
hereafter we call this approach Fundamental Plane 
limited (FPL). In the second recipe, we introduce a 
model for prolonged black hole gas accretion physi- 
cally motivated by galaxy dynamics. In this model, 
the subgrid physics (stellar and AGN feedback, as well 
as accretion disk microphysics) is bundled into one pa- 
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rameter for gas accretion efficiency tliat can be con- 
strained by future small-scale, fully general relativistic 
magncto-hydrodynamic AGN simulations. With this 
approach, gas accretion onto the black hole is based on 
the dynamics of the galaxy merger - fast at high rate 
for major mergers and prolonged at low rate for more 
minor mergers. This model does not use the black hole 
fundamental plane or the M - tr relation to limit black 
hole growth. During a minor merger, the accretion 
time scale is longer, but the accretion rate is lower, 
since less gas is driven toward the BH by a smaller 
global perturbation. For major mergers, the accretion 
rate increases, though in practice it never reaches the 
Eddington rate. We call this approach Prolonged Gas 
Accretion (PGA). In addition to gas accretion, black 
holes grow through mergers with other black holes too. 
In both recipes, we model SMBH growth from direct 
mergers including gravitational recoil. Since the recoil 
velocity is highly dependent on BH spin and orienta- 
tion, wc create ~ 20,000 merger trees to examine the 
combination of BH spin parameters that favors the M 
- a relation at redshift zero. 

The key difference between FPL and PGA mod- 
els is in the implementation of AGN feedback. In FPL, 
the black hole fundamental plane is used to calibrate 
AGN feedback. In PGA, we use AGN feedback to sup- 
presses the efficiency of gas accretion, while calibration 
parameters come from detailed small scale simulations 
of galaxy merger remnants. 

We describe our method in section 2 and intro- 
duce our black hole growth prescriptions. In section 
3, we present results, focusing on the evolution of our 
Sgr A* analogue in the FPL and PGA models. In Sec- 
tion 4, we present predictions for the local population 
of massive black holes, rogue black holes in the Milky 
Way, and ejected black holes into the nearby intcr- 
galactic medium. We discuss the implications of our 
results and future work in section 5. 



2 METHOD 

2.1 VL-2 Dark Matter Halo Merger Tree 

In this paper, we use publicly-available VL-2 evo- 
lutionary tracks (Diemand et al. 2007,2008). VL-2 
is the highest resolution cosmological N-body sim- 
ulation of Milky Way formation and evolution, and 
evolves a ACDM universe with WMAP3 parameters 
(nM=0.238, Oa =0.762, c78=0.74 and h=0.7). The high 
resolution region of VL-2 is embedded within a pe- 
riodic box with a comoving length of 40 Mpc. The 
evolutionary tracks consist of the time evolution from 
z=27.54 to z=0 of ~ 20,000 dark matter halos identi- 
fied at redshift z=4.56. Halos have masses as small as 
- \if Mq. 

In our hybrid method, we combine dark matter 
halo merger trees obtained from VL-2 evolutionary 
tracks with an analytical treatment of the physical 
processes that arise in the dynamics of galaxy and 



black hole mergers. Our N-body approach stops with 
the creation of the halo merger tree. We seed dark 
matter halos with Population HI black holes until z=5, 
following the work of Trenti & Stiavelli 2009 (Fig 1 in 
Trenti et al. 2009), and follow their merger history 
from redshift 27.54 to by constructing numerical 
merger trees interpolating between the snapshots at 
z< 15. To define the structure of each dark matter 
halo within the N-body generated merger tree, we as- 
sume a Navarro, Frenk, &: White (hereafter, NFW) 
density profile (Navarro, Frenk & White 1996). We 
set the parameters of a given NFW halo using the ap- 
proach presented in Bullock ct al. 2001, assuming the 
typical virial mass of a dark matter halo to be Mtyp 
= 1.5 X 10" M© at redshift zero. Note that only ha- 
los that are seeded with black holes or merging with 
other halos are retained in our merger tree. For each 
merger, we tag the more massive halo as the primary, 
with mass Mp, and the less massive halo as the sec- 
ondary or satellite halo with mass Ms. Note that the 
properties of low mass dark matter halos in the mass 
range of VL-2 (10^ Mq < M < 2x10^^ M©) have not 
been studied in detail at z > 3. 

Note that VL-2 evolutionary tracks do not include 
any halos or subhalos that were not present at z=4.56. 
Some halos are completely disrupted in mergers prior 
to that time (no subhalo counterpart at z=4.56), and 
some halos form after z=4.56. These halos arc small 
in mass and would have been cut out of the merger 
tree anyway since they are either too small to host the 
black hole seed or too small to have relevant amounts 
of cold gas for either black hole accretion or star for- 
mation. Note that VL-2 does not necessarily repre- 
sent the actual merger histories of Milky Way or Local 
Group. 

2.2 SMBH Growth Prescriptions 

We test two fundamentally different models for SMBH 
growth: Fundamental Plane Limited (FPL) and Pro- 
longed Gas Accretion (PGA). In FPL, gas accretion 
onto the black hole is controlled by AGN feedback in 
such way that the SMBH mass docs not grow above 
the black hole fundamental plane at any time. In PGA, 
we bundle star formation and AGN feedback, as well 
as accretion disk microphysics together into gas accre- 
tion efficiency parameter; the mass is not restricted to 
lie on the black hole fundamental plane. 

Our only guideline is the mass of Sgr A* observed 
today. Since the mass of the SMBH at the center 
of the Milky Way is well-constrained, we can iden- 
tify which physical processes during the growth are 
most important for the final black hole mass and test 
their parameter space. The first of these parameters 
applies to both models: the black hole seed's initial 
mass function (BHMF). We test various BHMFs and 
various values for the minimum and maximum black 
hole seed mass. A second effect that is present in both 
models is gravitational wave recoil (see Section 2.4). 
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The third parameter is unique in each model. In the 
FPL model, the third parameter is the mass accreted 
during the "radio" or quiescent mode of the AGN duty 
cycle. The SMBH growth during this phase is usually 
approximated by Bondi-Hoyle gas accretion, though 
this is not strictly necessary. We test various values for 
average mass accretion rate during this radio mode. 
Since the Bondi-Hoyle mechanism is commonly used, 
we may occasionally address this quiescent part of 
SMBH growth as the Bondi-Hoyle phase, but our ap- 
proach does not depend on the actual growth mech- 
anism here. In the PGA model, the third parameter 
is the gas accretion efficiency. In our implementation 
this parameter contains information on how the micro- 
physics acts to suppress accretion onto the black hole, 
such as SNe feedback, accretion disk physics, etc. We 
treat all of these processes as free parameters to study 
the observable consequences at z—Q. 

In all models, we adopt Gncdin (2000) and 
Kravtsov et al. (2004) initial cold gas fractions (fcg 
~ 10-5 for Mhaio ~ 3x10® Mo; fcg ~ 10"^ for Mhaio 
~ 4x10^ Mo; fcg ~ 0.13 for Mhaio > IxlO" Mq) in 
high redshift galaxies for a reionization range of 10 
< z < 11. In this manner we follow the approach of 
Somerville et al.(2008) for creating the gas reservoirs 
used for star formation and gas accretion onto the 
black holes. We do not include further multi mode 
accretion in the merger tree (i.e. Keres et al. 2005). 

2.2.1 Fundamental Plane Limited Model 

In the first model, FPL, we use the approach adopted 
by Somerville et al. (2008), but focus only on those as- 
pects that are important for SMBH growth. In FPL, 
the SMBH grows through mergers and gas accretion. 
Gas accretion is controlled by two modes of AGN feed- 
back. The quasar or "bright" mode occurs after the 
host galaxy goes through a major merger; this causes 
the central black hole to accrete gas exponentially un- 
til it reaches the Eddington limit, and it continues 
accreting at Eddington limit until it reaches a mass 
specified by the black hole fundamental plane. By de- 
sign, the black hole mass at the end of each merger is 
proportional to the stellar mass of the spheroid Mgph 
in the model: 

log(AfBH/Af,ph) = -3.27+0.36 erf [(/gas -0.4)/0.28], (1) 

where /gas is the amount of cold gas available for ac- 
cretion (Somerville et al. 2008). The quasar mode is 
replaced by a radio mode until the next major merger. 
During the radio mode, the black hole accretes gas 
at lower rates associated with radiatively inefficient 
Bondi-Hoyle accretion, such as what occurs in an ad- 
vection dominated accretion flow model (Fender et al. 
2004). The Bondi-Hoyle accretion rate controls the 
black hole mass during the longer quiescent phase be- 
tween galaxy mergers. In the Milky Way today, the 
SMBH mass is 4.2x lO'^Mo (Ghez et al. 2008), and the 
Bondi-Hoyle rate is ~ 5 x 1Q~^ Moyr"^ (Quataert & 



Gruzinov 2000, Melia & Falckc 2001). Since we know 
the current Sgr A* mass and the accretion rate ob- 
served today, we can determine how well the FPL 
model can match observations. To simplify our inves- 
tigation we adopt an average Bondi-Hoyle accretion 
rate during the radio mode over a Hubble time. The 
Bondi-Hoyle accretion rate in our model is averaged 
over all radio mode phases the galaxy goes through 
and since it explicitly excludes AGN feedback which 
would blow away the gas, the Bondi-Hoyle accretion 
rate we use is a lower constraint. If we model AGN 
feedback, the derived accretion rates would be higher. 

2.2.2 Prolonged Gas Accretion Model 

The SMBH in our PGA model grows through a com- 
bination of black hole mergers and gas accretion, in- 
spired by Micic et al. (2007). To review the approach, 

we assumed that major galaxy mergers would fun- 
nel gas to the black hole in each of the progenitors 
and activate an Eddington-limited growth phase for 
a Salpeter time. The black hole mass in Micic et al. 
(2007) grew as: MBn(t) = MBH.o(to) exp(At/tsai), 
where At = t - to, tsai = eMBHC^/[(l — e)L], e is 
the radiative efficiency, L is the luminosity, and c is 
the speed of light; in this picture the black hole mass 
would roughly double in 40 Myr (Hu et al. 2006). We 
distinguished two cases depending on the mass ratio of 
merging dark matter halos. The first is a more conser- 
vative criterion that allows black holes to accrete gas 
if the mass ratio of the host dark matter halo is less 
than 4:1 - a major merger. The second case sets an 
upper constraint on the final black hole mass by allow- 
ing gas accretion as long as the merging dark matter 
halos have a mass ratio less than 10:1 - on the cusp 
of what is considered a minor merger. Since our black 
holes merged promptly after the halos merged, the ac- 
cretion timescale and efficiency for major mergers was 
the same regardless of the mass ratio or redshift. 

In this paper, we continue to model the black 
hole growth as one of extended gas accretion excited 
by major mergers. At high redshift, this is likely a 
good assumption, though note that at low redshift 
when mergers are infrequent, secular evolution such as 
bar instabilities may dominate the gas (and therefore 
black hole) accretion. Integrated over the whole of a 
black hole lifetime, though, this major merger-driven 
gas accretion is likely to be the dominant source of 
gas inflow. Since the black hole growth is so strongly 
dependent on what fuel is driven to the center during 
galaxy mergers, it is important to characterize this 
merger-driven mass growth, including the critical gas 
physics that may inhibit or strengthen this nuclear 
supply. We are motivated by numerical simulations 
that include radiative gas cooling, star formation, and 
stellar feedback to study the starburst efliciency for 
unequal mass ratio galaxy mergers (Cox et al. 2008), 
which finds that the gas inflow depends strongly on the 
mass ratio of the galaxy (see also, e.g., Hernquist 1989, 
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Mihos & Hernquist 1994). This study parametrizes 
the efflciency of nuclear star formation (i.e. gas sup- 
ply and inflow), a, as a function of galaxy mass ratio. 
In a broader sense, this study shows how much gas 
is available for either star formation or gas accretion 
onto the central black hole. These two processes com- 
pete for the same gas and the outcome (nuclear clus- 
ter versus SMBH) depends on the mass of the host 
spheroid (Nayakshin et al. 2009). The efficiency of gas 
inflow is described by: 



a = aslope 



Ms 
Mo 



(2) 



where ao defines the mass ratio below which there is 
no enhancement of nuclear star formation (i.e. gas in- 
flow), and aslope is the fitted slope of the solid line 
in Cox ot al. (2008), Fig 15. Here, the gas accre- 
tion efficiency has a maximum of 0.56 for 1:1 halo 
mergers with asiopc=0.6, and falls to zero at ao. This 
parametrization is insensitive to the stellar feedback 
prescription. We use a to define how efliciently the 
merger funnels the galaxy's gas to the black hole ac- 
cretion disk. In general, all feedback processes (stellar 
or black hole) can be contained in agiope- Only small 
scale numerical simulations of gas accretion in AGNs, 
with fully implemented and resolved gas physics can 
relate Osiope to the mass ratio of merging galaxies. In 
Cox et al. (2008), the fitted ao parameter suggests gas 
infiow is sharply curtailed for mass ratios larger than 
9. We set the cut off mass ratio to Mp/Ms = 10 in 
order to compare with our previous black hole growth 
prescription (Micic et al. 2007). We adopt Qsiope = 0.6 
as in Cox et al. (2008). 

Now that we have implemented a realistic descrip- 
tion of the merger time for each black hole within a 
halo (see next section), we allow them to grow for a 
physically-motivated accretion timescale. The accre- 
tion of gas onto both the incoming and central black 
hole starts when the two black holes are still widely 
separated, at the moment of the first poriccntor pas- 
sage, and continues until the black holes merge (c.f. 
Di Matteo et al. 2005, Colpi et al. 2007). This sets the 
accretion timescale, tacc, as follows: tacc = tdf (r=rvir) - 
tdyn(r=rvir), wherc tdt(r=rvir) is the merger timescale 
including dynamical friction, and tdyn(r=rvir) is dy- 
namical time at virial radius rvir, which marks the 
first pericenter pass of the black hole. By stopping the 
accretion as the black holes merge, we roughly model 
the effect of black hole feedback in stopping further 
accretion. 

Putting these pieces together, the mass accreted 
by a black hole during tacc(r=rvir) is: 



Ma, 



MBH,o(e '-r - 1), 



(3) 



where Mbh.o is the initial black hole mass, a is the 
starburst efficiency (Cox et al. 2008), and tsai is de- 
fined above. After tat, the incoming black hole merges 
with the SMBH at the center and a new SMBH is 



formed after having accreted gas for tacc. The ac- 
cretion time and efficiency both implicitly encode 
the large-scale dynamics of the merger and the bulk 
gas accretion into the nuclear region, while tsai de- 
scribes the accretion disk physics. As before, wc set 
the Salpeter accretion efficiency to 0.1 (Shakura & 
Syunyaev 1973). 



2.3 Dynamical Priction 

Dynamical friction allows massive black hole binaries 
to form at the center of a galaxy in two ways. First, 
dynamical friction expedites the merger of two dark 
matter halos and later the merger of the galaxies they 
host. In this manner, merging galaxies can efficiently 
shepherd massive black holes to the center of the now 
system, roughly to the inner kiloparsec (see Colpi et 
al. 2007 for a review). Second, dynamical friction from 
the gas in the disk carries black holes deeper toward 
the galactic center, where they form binary and even- 
tually merge (e.g. Bcgclman, Blandford & Rccs 1980, 
Escala et al. 2005, Kazantzidis et al. 2005, Dotti et al. 
2007). We model both effects as foUows: 

The time for a satellite to sink to the center of 
a primary can be approximated using Chandrasekhar 
dynamical friction (Binney & Tremaine 1987): 



tc 



1.17rgi,,vc£" 
InA GMs ' 



(4) 



where InA is the Coulomb logarithm, InA ~ 
In (1 -|- Mp/Ms). To define the satellite orbit, we adopt 
values suggested by numerical investigations (Colpi et 
al. 1999) and used in previous semi-analytical work 
(Volonteri ct al. 2003): the circularity e° = 0.8, and 
the circular velocity Vc is determined at rdrc = 0.6 
rvir. For a 10^^ M@ halo, rvir ~ 300 kpc. 

Assuming that each merging galaxy carries a mas- 
sive black hole at its center, tfric is the merging time 
for massive black holes when all other processes (3- 
body scattering, gas dynamical friction, gravitational 
radiation, etc.) involved in the formation and later 
shrinking of the black hole binary are efficient and fast. 
Due to tidal stripping and possible resonant interac- 
tions, simulations have shown that the Chandrasekhar 
formula underestimates the merger time, especially in 
the case of minor mergers (e.g. Holley-Bockelmann & 
Richstone 1999, Weinberg 1989). If this is true, then 
semi-analytic studies of black hole merger rates using 
a Chandrasekhar formalism for the merger time will 
overestimate the true number of black hole mergers. 

In an effort to better parametrize dynamical fric- 
tion, Boylan-Kolchin ct al. 2008 used N-body sirrm- 
lations to study dark matter halo merging timescales, 
and confirmed that the Chandrasekhar formalism does 
underestimate the merger time, by a factor of fs 1.7 
for Mp/Ms ~ 10 and a factor of ~ 3.3 for Mp/Ms « 
100. They propose a fitting formula that accurately 
predicts the timescale for a satellite to sink from the 
virial radius to the host halo center: 
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Figure 1. Time for each satellite to reach the primary 
halo center as a function of the mass ratio of merging ha- 
los. Black circles represent the Chandrasekhar dynamical 
friction time and red pluses represent the merging time 
calculated from a simulation-based numerical fit (Boylan- 
Kolchin 2008). Both timescales are compared to the Hub- 
ble time and mass ratios are compared to major mergers. 
Merging halos above the horizontal black line will not fin- 
ish their merger and are removed from the merger tree, but 
their positions are updated within the primary halo. 



T"dyn 



A 



ln(l-FMhost/Msat) 



exp 



GA/ho 



1/2 



(6) 



where b = 1.3, c = 1.9, d = 1, A = 0.216, circularity j 
/ jc(E) = 0.5, rc(E) / rvir = 0.65, as defined in Boylan- 
Kolchin et al. 2008. The dynamical time, tdyn, is given 
at virial radius as: 



where Vc(rvir) = (GMvir/R)^''^ 

Figure 1 shows the dark matter halo merger 
timescale for each pair of merging halos for both dy- 
namical friction estimates. We have calculated the 
dark matter halo merger rate with both the Chan- 
drasekhar dynamical friction formula and the Boylan- 
Kolchin numerical fit (HoUey-Bockelmann et al. 2010). 
Unless otherwise stated, we adopt the Boylan-Kolchin 
timescale for all calculations. 

It is expected that in gas-rich galaxies, dynamical 
friction from the gas would bring two black holes close 
enough to form a binary whose orbit would shrink ef- 
ficiently, passing quickly from a binary in the 3-body 
scattering phase to one emitting significant gravita- 
tional radiation (Escala et al. 2005, Kazantzidis et al. 
2005, Dotti et al. 2007). Numerical simulations indi- 
cate that two black holes will sink from ~ 1 kpc to 
form a binary with a separation of less than a parsec 
in ~ 10 Myrs. We incorporate this physics by calcu- 
lating the dynamical friction timescale from the virial 



radius to the inner kpc, and then assume that the two 
black holes merge 10 Myr afterward. In practice, the 
dynamical friction timescale from equations 1 and 2 
from the inner kpc to the bound binary stage is often 
of order 10 Myr; the power of this gas-rich assumption 
lies in that it entirely circumvents the so-called 'final- 
parsec' problem thought to exist for low mass ratio 
mergers of lO^ '^ M© < Mbh < 10* M© within static, 
spherical, gas-poor galaxy models (e.g. Milosavljevic 
& Merritt 2003). We explicitly assume that the black 
holes in our simulation do not stall at the final parsec 
before merger, and instead are ushered efficiently into 
the gravitational radiation stage where they coalesce; 
this assumption may hold even in the case of gas-poor 
galaxy models as long as the model is not spherical or 
in equilibrium (e.g. HoUey-Bockelmann & Sigurdsson 
2006, Berczik et al. 2006, R. Spurzem, private com- 
munication). Note that this implies that if we assume 
that each halo initially carries a black hole at its cen- 
ter, and that each host galaxy is gas-rich, figure 1 also 
estimates the time for a black hole to sink from the 
virial radius to the center of the host galaxy, become 
a bound black hole binary, and inspiral due to gravi- 
tational radiation. 

In our initial work (Micic et al. 2007), mergers 
of dark matter halos trigger the immediate merger of 
the black holes they are hosting. In this paper, sub- 
sequent mergers of the central black holes are delayed 
to account for dynamical friction of the halos and the 
black holes within the galaxy. Black holes will not 
, (5^erge if their merger time is larger than a Hubble 
time, and in that case, we advance the black hole posi- 
tion within the primary halo at each timestep. Know- 
ing the dynamical friction timescale for each merger, 
we postpone the black hole mergers accordingly. For 
the final kpc, we assume that ambient gas and/or non- 
sphericity will cause two black holes to coalesce within 
10 Myrs. 



2.4 Gravitational Recoil 

Binary black holes strongly radiate linear momentum 
in the form of gravitational waves during the plunge 
phase of the inspiral - resulting in a "kick" to the new 
black hole. This, in itself, has long been predicted as 
a consequence of an asymmetry in the binary orbit or 
spin configuration. Previous kick velocity estimates, 
though, were either highly uncertain or suggested that 
the resulting gravitational wave recoil velocity was rel- 
atively small, astrophysically speaking. Now, recent 
results indicate the recoil can drive a gravitational 
wave kick velocity as fast as ~4000kms~^ (Herrmann 
et al. 2007, Baker et al. 2007, Gonzalez et al. 2007a, 
2007b, Koppitz et al. 2007, Campanelh et al. 2007, 
Schnittman & Buonanno 2007, Lehner & Moreschi 
2007, McWilliams 2008, Lousto & Zlochower 2009, 
Sperhake 2009). In reality, much smaller values than 
this maximum may be expected in gas-rich galaxies 
due to the alignment of the orbital angular momentum 
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and the spins of both black holes (Bogdanovic et al. 
2007). Q 

Recent studies also hint at a potential purely 
general relativistic spin alignment mechanism (Kes- 
den et al. 2010). However, even typical kick velocities 
(~ 200 kms^^) are interestingly large when compared 
to the escape velocity of most astronomical systems - 
low mass galaxies, as an example, have an escape ve- 
locity of ~ 200 kms~^ (e.g. Holley-Bockelmann et al. 
2007). The effect of large kicks combined with a low 
escape velocity from the centers of small dark matter 
halos at high redshift may play a major role in sup- 
pressing the growth of black hole seeds into SMBHs. 
Even the most massive dark matter halo at z>ll can 
not retain a black hole that receives > 150 kms~^ kick 
(Merritt et al. 2004, Micic et al. 2006). We incorpo- 
rate the effect of recoil velocity on the growth of the 
SMBH by assigning a kick to each black hole merger 
in our merger tree. We follow the approach adopted 
by Holley-Bockelmann et al. (2007), which uses the 
parametrized fit of Campanelli et al. 2007 to general- 
ize the recoil velocity as a function of the mass ratio of 
the merging black holes, each individual black hole's 
spin amplitude, and the alignment to the orbital an- 
gular momentum. We assume that as black holes form 
a hard binary, their orbits will be highly circular and 
the eccentricity is close to zero: 

2ll/2 



where 



= A 



1) 



and 



(1 + qf 

if cos (e - eo) 



i + s- 



(7) 



(8) 



(9) 



+ 



Q2 



(10) 



Here, the fitting constants are A= 1.2 x 10* kms 



(7.3 ± 0.3) X lO'' kms" 



and 



B= -0.93, H 

Kcos(e - eo) = (6.0 ± 0.1) X 10* , while the sub- 
scripts 1 and 2 refer to the first and second BH re- 
spectively; _L and || stands for perpendicular and par- 
allel to the orbital angular momentum; the mass ra- 
tio q= M2/M1; the reduced spin parameter Ui = 
Si/Mf, where Si is the spin angular momentum of 
BH i. The orientation of the merger is specified by: 
G, the angle between the "in-plane" component of 
5' = (Ml + M2) (SI/M2 - SI /Ml) and the infaU di- 
rection at merger; Go, the angle between 5' and the 
initial direction of motion; and the angle between 
the unequal mass and spin contribution to the recoil 
in the orbital plane. 

We assume that the orbit is circular, and we take 
the mass ratio of merging black holes directly from our 



^ However, recent observations (Komossa et al. 2008) may 
sliow evidence for large kick velocities ~ 2500 kms^^. 
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Figure 2. Distribution of gravitational recoil as a func- 
tion of the mass ratio of merging black holes on circular 
orbits. The black lines represent random spin orientations 
and amplitudes - the Kl distribution. In red is the kick 
range when the spin orientation of the binary black holes 
aligns with the orbital momentum of the binary - the K2 
distribution. If the kick is larger than the escape velocity, 
the black hole is ejected from the halo. 



merger tree. We have the following free parameters: 
the spin amplitude and orientation of each black hole, 
and the orientation of the merger. We explore two spin 
distributions. The Kl model chooses the spin param- 
eters from a uniform distribution, while the K2 model 
assumes the black hole spins are aligned with the or- 
bital angular momentum (Bogdanovic et al. 2007). 
Figure 2 shows the distribution of kick velocities for 
both models. We apply both kick velocity distribu- 
tions to the mergers in our merger tree using 1000 
realizations each. This yields 2,000 Milky Way merger 
trees for the PGA model and 2,000 for the FPL model. 

For each black hole merger, we calculate the kick 
velocity and compare it to the escape velocity of the 
host halo at the time of merger. If the kick is larger 
than the escape velocity, the resulting black hole is re- 
moved from the center of the host halo and from the 
merger tree. As the result, the host halo will not have 
a central black hole for a period of time. Subsequently, 
when the next black hole sinks to the center, it will 
simply take place of the central black hole. Depend- 
ing on how often the central black hole is ejected, the 
final SMBH mass may be substantially smaller. The 
sequence and number of kicks in one merger tree real- 
ization depends on our kick distribution. Over many 
realizations, we can produce a probability function for 
the final black hole mass for each model. 
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Figure 3. Thick black - the average change in the merger 
redshift when dynamical friction is applied as a function of 
merger redshift (excluding dynamical friction). Red lines 
represent the minimum and maximum shift in the black 
hole merger redshift. Dynamical friction postpones black 
hole mergers toward lower redshifts, making them louder 
LISA sources and increasing the local rate as well. 

3 RESULTS 

Black hole mergers postponed for longer than a Hub- 
ble time will not occur, which reduces the merger rate 
over all redshifts. Those mergers that are postponed 
but have a merging timescale less than a Hubble time 
will merge at lower redshifts. This is made explicit in 
figure 3, which plots the average change in redshift for 
a merger as a function of redshift if one includes dy- 
namical friction. Black hole mergers at redshift 10 are, 
on average, pushed to redshift 7, for example. This re- 
sults in an increase in the black hole merger rate at low 
redshifts. Since each merger occurs at a lower redshift, 
it will certainly be a louder gravitational wave source, 
and any associated electromagnetic signature will be 
brighter, as well (HoUey-Bockelmann et al. 2010). 

Interestingly, including dynamical friction from 
gas at distances smaller than 1 kpc does not make dif- 
ference in our merger time scales. Those black holes 
that reach 1 kpc from the galactic center will reach a 
parsec even without the gas - gas is simply invoked 
to ensure that the black holes pass through the 3- 
body scattering stage efficiently. For the rest of this 
paper, we adopt the Boylan-Kolchin dynamical fric- 
tion fit as a realistic treatment of the black hole merger 
timescale. 

3.1 Realizations of Supermassive Black Hole 
Growth 

Figure 4 shows the merger rates in our volume for 
those dark matter halos that merge in less than a Hub- 
ble time within the mass range 10^ M0 < Mdmh < 
10^^ Mq, for various mass ratios and combined masses 
of merging dark matter halos. In our volume, all the 
mergers that occur at redshifts z < 1.5 will not fin- 
ish in less than a Hubble time, as these mergers axe 




Figure 4. Dark matter halo merger rates as a function of 
redshift for all halos that finish merging by redshift zero. 
Three panels show various ranges for the halo mass ratio 
and the total combined halo mass. Mp is the primary halo 
mass, Ms is the satellite mass, and the total mass is defined 
by p = log (Mp -I- Mg). The panel with Mp/Mg < 10 
shows merger rates for major mergers that activate gas 
accretion onto their black holes. 



all high mass ratio. We remind the reader, though, 
that these halo merger rates are calculated for VL-2 
40 Mpc'^ simulation box. This model, by design, has 
no Galaxy Clusters, or even rich local groups, and it is 
for sparse local groups. The VL2 host halo was specif- 
ically chosen "not" to have undergone a major merger 
more recently than about z~l. Almost all the low red- 
shift mergers for this mass range SMBH are in richer, 
denser environments, which are just not covered in 
this model. A better global rate for these low mass 
halos will be achieved with more of these volumes. 

Since in PGA the merger timescale is larger for 
higher mass ratio halos, the black hole gas accretion 
time scale is longer as well. Figure 5 shows this time 
scale for black holes hosted by dark matter halos that 
merge with mass ratio Mp/Mg and combined mass 
p=log(Mp-|-Ms). Notice that gas accretion is activated 
only for a small range of halo mass ratios Mp/Ms < 
10 (the dashed line in figure 5). For these mergers, the 
black holes do not accrete for longer than ~ 1000 Myr, 
and the accretion efficiency is strongly damped as the 
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Figure 5. Gas accretion time, tacc, for a central black hole in a satellite halo of mass Msat that is merging with a primary 
halo of mass Mp^inj in the PGA model. The total mass is denoted by p = log (Ms +Mprim), and tacc = tmerger - tdyn where 
tmorgcr is the merger time scale, and t^y^ is the dynamical time, making accretion begin at the first pericenter pass. Halos 
that merge with mass ratios less than 10 (vertical line) will allow gas accretion onto the black holes. Only satellite halos 
with mass < lO^'^ Mq will have an accreting central black hole. 



halo mass ratio tends toward 10:1 - in no case does 
the system reach the Eddington hmit. 

Figure 5 oflters valuable information about the 
FPL model as well. Panels e and f show that there are 
no mergers with mass ratios 10:1 or less for halos above 
lO^'^ Mq. In fact, the halo that VL-2 marks as the 
Milky Way halo does not go through any major merg- 



ers at low redshift. In the context of FPL, the quasar 
mode of rapid gas accretion occurs in Milky Way pro- 
genitors at high redshifts while the radio mode dom- 
inates most of the Sgr A* growth at low redshifts. 
We should also point out that VL-2 merger history of 
Milky Way is just one of the possible merger histories, 
and perhaps not the true one. 
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3.1.1 The Initial Mass Function of Black Hole 
Seeds. 

Galaxies in the mass range of the VL-2 simulation 
have values for escape velocities comparable to the 
amplitudes of gravitational wave recoil for black hole 
mergers with mass ratios q>0.1, even in the case of 
the lower K2 kick distribution. If we are to seed these 
galaxies with their first black holes, the choice of the 
black hole mass function (BHMF) will determine the 
mass ratio of merging black holes. We examine the in- 
fluence of the choice of BHMF on VL2 merger trees 
in 1,000 kick realizations for BHMFs with constant 
values of 100; 1000; and 5000 Mq seeds; and also for 
random values in the mass ranges 10 - 200 Mq; 10 
- 500 Mq; and 10 - 1,000 Mq. We find that having 
a constant BHMF or one in a narrow mass interval 
leads to close-to-equal mass ratio black hole mergers. 
These mergers have the highest gravitational wave re- 
coil and ejects of black holes in high redshift dwarf 
galaxies, which suppresses SMBH growth. Due to the 
uncertainty in the formation channel for seed black 
holes, as well as the lack of theoretical constraints on 
primordial gas fragmentation during PopHI star for- 
mation (Glover et al. 2010), 10 - 1,000 Mq is the least 
biased of our BHMFs. We adopt this flat BHMF here- 
after for the seeds in our merger trees in both the FPL 
and PGA models. Note that the seed mass in the FPL 
model is much less relevant than in the PGA model, 
because the black hole can not grow more massive 
than what is dictated by the black hole fundamental 
plane, regardless of the initial seed mass. 

3.1.2 Massive Black Hole Accretion Properties 

With the adopted flat BHMF, the only constraints in 
the FPL model are the time-average Bondi-Hoyle ac- 
cretion rate during the radio mode and the strength of 
gravitational wave recoil. We start by running Milky 
Way merger trees for various values of the accre- 
tion rate and recoil distribution, looking for the val- 
ues which, in combination with the accretion during 
quasar mode, produce a ~ 4.2 xlO'' M© black hole at 
z=0 in the majority of merger tree realizations. This 
condition is satisfled for the merger trees with rribondi 
= 7 xlO~* M0yr~^ (brown horizontal line in figure 
6 up), consistent with the observations of accretion 
rates in low mass AGNs in the local Universe. It has 
been suggested that the Milky Way has evolved from 
a similar low mass AGN, which explains the fact that 
the observed value for accretion in Sgr A* today is an 
order of magnitude smaller than this best fit accretion 
rate. Outside of the Sgr A* analogue, there are also 
35 field dwarf gala:xies at z=0 with central black holes 
(see section 4.1). Figure 6 shows black hole accretion 
properties during the quasar and radio mode for the 
merger tree that excludes recoil, which corresponds 
to the maximum accretion. In this model. Milky Way 
progenitors go through their last quasar mode between 
z=2 and z—3 (upper panel. Fig 6), after which Bondi- 




Figure 6. Accretion rates and luminosities during the 
quasar and radio modes for the FPL model (recoil ex- 
cluded). The upper panel shows gas accretion rates and 
luminosities for all black holes during a random short time 
interval {St = tf - ti). Points represent the quasar mode, 
while the brown horizontal line represents the radio mode. 
The black hole mass in the quasar mode increases at lower 
redshifts, with the final accretion episode between z=3 and 
z=2. Bondi - Hoyle accretion is the dominant growth mech- 
anism at redshift z < 2. The bottom panel shows the gas 
accretion rate at short time intervals compared to the Ed- 
dington rate (quasar and radio modes combined). Most 
black holes accrete at the Eddington rate by design. 



Hoyle accretion becomes the dominant mechanism of 
Sgr A* growth (horizontal line in upper panel, Fig 6). 
In the bottom panel of figure 6, the same accretion 
rate is compared to the Eddington rate. Most of the 
accretion occurs at the Eddington rate - this is ex- 
pected since the assumption is built in that the black 
holes reach the Eddington rate very quickly. In many 
cases, however, accretion is not activated at all be- 
cause the host galaxy is small and the spheroid mass 
prevents the growth of central black hole to match the 
black hole fundamental plane. 

By excluding kicks, figure 7 shows the upper con- 
straint on black hole accretion properties for all black 
holes in the PGA merger tree. Recall that the BHMF 
is fiat and the accretion efficiency is set by simulations. 
With the choice of BHMF and gas infiow efficiency 
fixed in this model, kicks are the only free parameter. 
Black holes in the PGA model can grow beyond the 
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Figure 7. Accretion in the PGA model. Left panels show black holes with mass less than 10* Mq, and right panels show 
more massive black holes. Upper panels: Gas accretion rates and luminosities for all accreting black holes during a random 
short time interval (St = tj - ti). Bottom panels: Gas accretion rates at short time intervals as a fraction of the Eddington 
rate for all accreting black holes. Unlike the FPL model, here most black holes accrete well below the Eddington rate. 
There are more (< 10* Mq) accreting black holes in the PGA model because their growth is not limited by the black hole 
fundamental plane. 



limit specified by the black hole fundamental plane. 
Instead of invoking AGN feedback in the form of an 
upper limit on the black hole mass, we incorporate it 
into the gas accretion efficiency. These feedback mech- 
anisms suppress the accretion efficiency, leading to ac- 
cretion rates that are 0.01 - 0.8 of Eddington (bottom 
panel. Fig 7). these lower accretion efficiencies also re- 
sult in lower luminosities in the PGA model. While 
for ~10'' Mq black holes, the maximum luminosities 
is ~ 10*^ L0 in the FPL model (upper panel. Fig 6), 
these black holes have an order of magnitude lower 
luminosity in the PGA model (upper panel. Fig 7). 
The maximum accretion rate of the PGA black holes 
at the redshift of their last accretion is ~ 0.01 - 0.1 
Mo/year or 0.8 rriEdd. Thereafter, the gas accretion is 
extremely damped. Although these SMBHs are far too 
low a mass to be considered a high luminosity AGNs 



analogue, it may be appropriate to link these SMBHs 
with low luminosity AGNs (Ho 2008). Our approach is 
in broad agreement with the observations of AGN life- 
times (Hopkins & Hernquist 2009), which show that 
most local black holes must have gained their mass in 
no more than a couple of accretion episodes. Again, 
we point out that Hopkins & Hernquist 2009 discuss 
their results for more massive black holes (>1O^-^M0) 
than those in our study. 



3.1.3 Black Hole Growth with Gravitational Wave 
Recoil 

We include gravitational wave recoil from the uni- 
form Kl and lower recoil K2 distributions and repro- 
duce 1,000 merger trees for each kick distribution in 
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Figure 8. The evolution of Sgr A* in tlie FPL model 
for 1000 Kl kick realizations. The black solid line shows 
the black hole growth when kicks are excluded. Solid and 
dashed red lines show the black hole growth that pro- 
duces the largest (4.3xlO^M0) and smallest (2.0x10'5Mo) 
final black hole. The most common black hole mass is 
3.9x10^ M0 (blue). Also in blue: the l-cr spread around 
the most common value for redshifts z=1.5; 1.0; 0.5; and 
0.0, further showing that the most common mass is near 
the mass of Sgr A*. 



both the FPL and PGA models. This results in 4,000 
merger trees presented in figures 8, 9, 10 ,and 11. 

Figure 8 shows results for the Kl distribution in 
the FPL model. When kicks are excluded, the final 
black hole mass is 5 x 10® Mq. The largest and small- 
est black holes at the center of the Milky Way analog 
are 4.3x10*^ M© and 2.0x10*^ Mq, respectively, while 
the most common black hole mass is 3.9 x 10® Mq. Fig- 
ure 8 also shows the 1-a spread around the most com- 
mon value for redshifts z—1.5; 1.0; 0.5; and 0.0, further 
showing that the typical mass remains near the mass 
of Sgr A*. In the realization that favors the smallest 
SMBH mass, the kicks are so large that the black hole 
is ejected from the center of the Milky Way in two 
occasions. The reason the FPL model produces a final 
SMBH mass that so closely matches Sgr A* is quite 
trivial: it is built into the model that any major merger 
will grow the SMBH up to the mass dictated by the 
black hole fundamental plane. Even when the central 
black hole in the Milky Way analog is only ~ IOOMq 
at z=l.l, there is plenty of time to grow the black hole 
back to the supermassive range. 

This also acts to shift the typical black hole mass 
toward the merger tree realization with largest final 
black hole mass. In the case of lower amplitude K2 
distribution (with black hole spin axes aligned) this 
effect is even more pronounced (figure 9). In figure 9, 
we show the z < 2 part of the merger tree for the 
case that excludes kicks, as well as the extrema in 
final black hole masses, all of which are close to the 
observed Sgr A* mass. 

Gravitational wave recoil matters much more in 




red.shift 

Figure 9. Sgr A* evolution in the FPL model for 1000 K2 
kick realizations (lower recoil). Black : the z < 2 part of 
the merger tree when kicks are excluded. Solid red: kick 
realization favoring the largest black hole at the center of 
Milky Way analog. Dashed red: the kick realization favor- 
ing the smallest final black hole. Both the smallest and 
largest values are close to the observed Sgr A* mass. 
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Figure 10. Sgr A* evolution in the PGA model for 1000 
Kl and K2 (lower recoil) kick realizations. Black solid line: 
Sgr A* when kicks are excluded. The red thick and thin 
lines show the Kl and K2 realizations favoring the largest 
black holes. The blue thick and thin lines show the Kl and 
K2 realizations favoring the smallest black holes. Here, the 
scatter is much larger than in the FPL model, particularly 
for the higher Kl recoil distribution. 



the PGA model, as can be seen in figure 10 which 
shows the central SMBH mass of the Milky Way ana- 
logue as a function of redshift. Although the final 
black hole mass lies in a much wider range of masses, 
these outliers in mass are rare. Figure 11 shows the 
black hole mass with the largest probability at red- 
shifts z=3.0 (blue); 1.5(red); and O.O(black) for PGA 
models. The black line shows that for the lower re- 
coil velocities, K2 case at z=0, the typical black hole 
mass is in the Sgr A* range 90% of the time. Note 
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Figure 11. Upper panel: Number of lower recoil K2 real- 
izations in the PGA model where Sgr A* reaches a certain 
mass at a certain redshift. Bottom panel: Number of Kl 
realizations in the PGA model where Sgr A* reaches a cer- 
tain mass at a certain redshift. Histograms are at z=3.0 
(blue dashed); z=1.5 (red thin); and z=0 (thick black). In 
the lower recoil K2 case, Sgr A* reaches the observed value 
at z=0 in more than 90% of kick realizations. 



that since it is thought that the black hole spin vec- 
tors may ahgn with one another before coalescence 
(Bogdanovic et al. 2007, Sperhake 2009, Kesden et al. 
2010), we consider this our best and most realistic 
model. If SMBH grows according to the PGA model 
then kicks must be in the lower range described by K2 
distribution. However, this does imply a scatter in the 
M -a relation at the low mass end. 



the satellite does reach the center of Milky Way, the 
following black hole merger might lead to the black 
hole ejection). 



4.1 Massive Black Holes in the Local Dwarfs 

For dark matter halo masses above roughly 5 x 
10^^ M0, the dark matter halo mass correlates well 
with the mass of the supermassive black hole (Fer- 
rarese 2002): 



Mbh 
108 Mp 



A\ 



Mdmh 
1012 M 







(11) 



where Mdmh is the dark matter halo mass, A=0.1 
, and m=1.65 defines the slope in the relation. Be- 
low this mass, there is tentative evidence that halos 
are less effective at forming massive black holes, and 
may even be unable to form them (Ferrarese 2002). 
To be fair, though, if this relation does hold down to 
a lO" M0 halo, the expected 2 x 10* M© black hole 
would be difficult to detect, observationally. 

Outside of the Milky Way, there are 35 local field 
dwarf galaxies in VL-2 identified at z=0 in the mass 
range lO'^ M© - lO" M©. Figure 12 shows the black 
hole mass at the centers of these field dwarfs as a func- 
tion of dark matter halo mass in the FPL (left) and 
PGA (right) models. Thick black vertical lines repre- 
sent minimum and maximum values for the black hole 
mass in 1,000 kick realizations. The Red pluses cor- 
responding to the vertical black lines represent most 
common mass between minimum and maximum. If a 
red plus is absent, the most common black hole mass 
for that halo is zero. Linear fits to the most common 
masses are presented by red lines. Note that for our 
linear fits, we exclude those field dwarfs that have 
ejected their central black holes entirely in more than 
50% of kick realizations. 

The FPL model predicts an A/bh-A/dmh for field 
dwarf galaxies in Kl and lower recoil velocities K2 
cases as follows. 

For FPL Kl: 
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(12) 



4 OTHER POPULATIONS OF MASSIVE 
BLACK HOLES 

We distinguish three different massive black hole pop- 
ulations in our merger trees: black holes at the centers 
of local field dwarf galaxies (other halos in VL-2 out- 
side of the Milky Way halo); rogue black holes scat- 
tered through the Milky Way halo (remnants of satel- 
lites which merge with Milky Way but do not reach 
the center in less than Hubble time, Micic et al. 2007, 
Bellovary et al. 2010); and black holes ejected from 
their host halos to occupy the intergalactic medium (if 



while FPL lower recoil K2 has a slightly shallower 
slope: 



1.06±0.41 
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(13) 



The PGA Kl model, on the other hand, predicts: 

/ ^ 0.64±0.36 

Mbh AfoMH 



108 M, 



© 



1012 M© 



(14) 



and finally, the PGA lower recoil K2 model has the 
following relation: 
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Figure 12. Black hole demography for the FPL Kl model - upper left; FPL lower recoil K2 model - bottom left; PGA Kl - 
upper right; and PGA lower recoil K2 - bottom right. The dashed line in every panel represents the Ferrarese relation with 
a slope = 1.65. Red pluses: the most common central black hole mass. For the rest of 35 field dwarf galaxies most common 
outcome is absence of central massive black hole. Red line: linear fit to the most common black hole masses. Black vertical 
bars: the scatter in final black hole mass produced by kicks. This would reflect the scatter in relations that connect the 
central black hole mass with the properties of field dwarf galaxies. 
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between FPL and PGA models might be insignificant, 
(15) considering tiiat scatter around the linear fits is sub- 
stantial. 



A generic prediction of each model is that the 
slope is shallower than the Ferrarese relation presented 
by dashed lines in Fig 12. For FPL Kl and lower re- 
coil K2 models, the slope m is 1.14 and 1.06, close to 
MuMU I AfBH = 10^. PGA models predict an even 
shallower slope of 0.64 in Kl and 0.80 in the case of 
lower recoil velocities K2. This tests the PGA model 
by predicting a substantially flatter mass function for 
lower mass SMBH. We do point that these differences 



4.2 Rogue Massive Black Holes in the Milky 
Way 

Rogue black holes are carried into the Milky Way halo 
by their host satellites, but have halo merger time 
scales longer than a Hubble time. Over the span of 
the simulation, the Milky Way accretes 669 satellite 
halos which host a black hole, and 249 have orbital 
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Figure 13. Number of rogue black holes spread through the Milky Way halo as a function of mass, for the FPL (left panel) 
and PGA (right panel) models. Red histograms - the merger tree that excludes black hole kicks. Blue histogram - the most 
common merger tree in the Kl kick distribution. Green histogram - the most common merger tree for the lower recoil K2 
case. These black holes originate from the centers of satellites that merged with the Milky Way analogue but had merger 
timescales larger than a Hubble time. As the black holes fail to reach the Milky Way center, they orbit in the Milky Way 
halo. These models predict hundreds of rogue black holes with masses reaching up to ~ 200,000 Mq. 



decay times that are longer than a Hubble time (see 
figure 1 and figure 4 for the merger time scale, merger 
rate and mass ratios of these halos). By redshift zero, 
then, these halos are still orbiting at significant dis- 
tances from the center of Milky Way - and although 
the halos themselves have been stripped by the pri- 
mary potential, they still host massive black holes. 
Figure 13 shows the number of rogue black holes as a 
function of mass for the FPL (left) and PGA (right) 
models. As before, the upper limit for black hole mass 
excludes kicks and is presented in red in both panels. 
For each kick model, we also present the number of 
rogue black holes with the typical mass over the real- 
izations. The PGA model predicts that there should 
be at least one massive rogue black hole in the Milky 
way halo with mass in range 10^ Mq -10^ Mq. 

The distance from the Milky Way center in all 
models is correlated with the redshift at which galax- 
ies merge (figure 14). If the merger occurred early at 
high redshift, the black hole had more time to sink 
to the center, and for mergers close to z=0, the black 
hole has just entered Milky Way halo. The redshift at 
which the halo merges with the Milky Way also marks 
the end of the rogue black hole growth. Hence, rogue 
black holes that enter the Milky Way at later times 
have had more time to grow, and will be more massive 
than those accreted through early mergers. Figures 
15 and 16 show how rogue black hole mass, satellite 
merger redshift, and distance from the Milky Way are 
correlated. Figure 15 represents the FPL model and 
figure 16 represents the PGA model. Kl realizations 
are shown in the left panels, and lower recoil K2 real- 
izations in the right. 

Unfortunately, both figures clearly show that the 
most massive rogue black holes are furthest away and 




r / R 



Figure 14. Redshift of satellite mergers with the Milky 
Way analogue as a function of the current position in the 
Milky Way halo. Those satellites that merged with the 
Milky Way at high redshift have sunk closer to the Galac- 
tic center. In later mergers, the satellites have just entered 
the Milky Way and are closer to the virial radius. Since 
each satellite potentially carries a black hole, this figure 
shows the distance of rogue black holes from the center of 
the Milky Way. 



may be very difficult to observe. We estimate the 
Bondi accretion onto these black holes to determine 
how observable they may be and we describe the pro- 
cess below. 

There is sufficient evidence that the halos of 
galaxies, from massive ellipticals to isolated spirals, 
can host diffuse hot halo gas (Matthews & Brighenti 
2003, O'SuUivan et al. 2001, Pedersen et al. 2006). 
As these black holes orbit within the primary halo, 
they can accrete from this ambient hot halo gas via 
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Figure 15. Rogue black hole mass as a function of Galactic distance and progenitor halo redshift for the FPL model. Left 
panels; the Kl distribution. Right panels: the lower recoil K2 case. The blue squares represent the most massive rogue black 
holes; red circles represent the typical black hole mass, and black dots represent the lightest. Those black holes carried by 
satellites that merge later have had more time to grow, so the rogue black hole mass correlates with both distance from the 
center and merger redshift. The most massive rogues reach ~ 30,000 Mq at 0.8 — 1.0 of virial radius. 



Bondi-Hoyle accretion. Here, the mass accretion rate 
can then be described by: 



(c2+l;2)3/2 ' 



(16) 



where pisM is the density of the halo gas, Cs is the 
gas sound speed, and v is the velocity of the black 



hole. For these rogue black holes, we add this more 
quiescent form of accretion from the time the black 
hole has entered the virial radius of the halo in or- 
der to determine whether any may be visible today. 
To model the hot halo gas, we use a technique out- 
lines in Sinha & HoUey-Bockelmann (2009). In short, 
we assume it has an isothermal density profile consis- 
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Figure 16. Rogue black hole mass as a function of Galactic distance and progenitor halo redshift for the PGA model. Left 
panels: the Kl distribution. Right panels: the lower recoil K2 case. The blue squares represent the most massive rogue black 
holes; red circles represent the typical black hole mass, and black dots represent the lightest. The most common rogue black 
hole mass remains at the initial black hole seed mass, but a number of massive rogue black holes are greater than 10,000 
Mq. The most massive is Mrogue ~ 230, 000 Mq at 225 kpc from the Galactic Center (0.8 - 0.9 of virial radius). 



tent with X-ray observations of halo gas in eUipticals 
(Matthews & Brighenti 2003), and that the ideal gas is 
in hydrostatic equilibrium with the gravitational po- 
tential of the halo. We assume a gas core radius of 
0.3Rs, where Rs is the scale radius of the NFW halo, 
and that the gas temperature is the virial temperature 
at the virial radius. This yields a central gas temper- 



ature of ~ 2 X 10*' if for a dark matter halo of mass 
2 X 10^^ Mq, which is also consistent with observations 
(Pedersen et al. 2006) . At redshift zero, we assume the 
black holes are on bound elliptical orbits of eccentric- 
ity 0.8 to match simulation predictions for satellite 
mergers (e.g. Ghigna et al. 1998, Sales et al. 2007). 

By redshift zero, these black holes will have grown 
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Figure 17. Bolometric luminosity of rogue black holes as a function of distance from the primary halo center for the FPL 
(left) and PGA (right) models. The Kl case is on the upper panels and the lower recoil K2 case on the bottom. Blue squares 
represent the most massive rogue black holes; red circles represent the typical black hole mass, and black dots represent 
the smallest. Black holes are assumed to accrete via a Bondi-Hoyle mechanism from the ambient gas in the Milky Way. 
Most rogue black holes are well below solar luminosity, but a few are expected to be luminous Xray sources (just under the 
ULX cut-off) that reside in the outskirts of the Milky Way halo. The most massive of them is in PGA model and has a 
luminosity of 4000 Lq . 

to a new mass, Mroguc, and will have a new accretion around a spinning black hole sets ~ 0.1 (Shakura 
rate from equation [161 The radiated luminosity from & Syunyaev 1973). We find that our most massive 
this process can be determined from LroguG = ??Mbhc^, rogue black hole Mroguc ~ 227,643 Mq, residing at 
where rj describes how efficiently mass is converted to 225 kpc from the Galactic center, can radiate at 3,834 
energy; the standard assumption for an accretion disk L0 (figure 17, PGA panels on the right), making it 
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Figure 18. Upper panels: The cumulative number of ejected black holes as a function of redshift for all halos in VL-2 for 
the FPL (left) and PGA (right) models. The Kl (solid line) and lower recoil K2 (dashed line) kick distributions are shown. 
In blue is the maximum cumulative number of ejected black holes; in red is the typical number; in black - the minimum. 
Lower panels: The black hole mass of the ejected population. The colors are the same as in the top panels. 



somewhat fainter than the brightest ULXs; the fact 
that the hot halo gas is so tenuous at these radii is 
what accounts for the relatively low luminosity. 

However, if these rogue black holes are common, 
they will be observable in the outskirts of the Milky 
Way halo with Chandra, and the most luminous can 
have sufficient signal-to-nose for spectral characteriza- 
tion. Note, though, that a significant fraction of these 



may have advection-dominated accretion flows, rather 
than thin accretion disks, and hence rj may be much 
less than 0.1; this would naturally make the black hole 
less luminous, but with a harder spectral signature. It 
has also been suggested that rogue black holes in the 
molecular or cold neutral part of the ISM should be 
more easily detected in the radio than in X-rays (Mac- 
carone 2004, Maccarone et al. 2005). Maccarone 2005, 
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for example, calculates that a 2,600 Mq black hole 
should be observable with LOFAR (above 30^Jy) in 
the radio out to about 15 kpc - much closer than the 
most massive rogue black hole in our model. Assum- 
ing larger black hole velocities and a more efficient 
Bondi accretion rate, rogue black holes may be ob- 
served at even larger distances. Even if they are not 
observable electromagnetically, they may be detected 
by mesolensing (di Stefano 2007). 



4.3 Ejected Massive Black Holes 

As the SMBH grows, incoming black holes are kicked 
from the center of the host halo, forming a popula- 
tion of "ejected" black holes. Figure 18 shows the cu- 
mulative number of black holes lost from the merger 
tree due to kicks larger than the escape velocity in 
the top panels, and the ejected black hole mass as a 
function of redshift on the bottom panels. Again, we 
show the FPL model on the left and our PGA model 
on the right. The redshift of the last SMBH ejection 
is between z=2-3 for the PGA model, and between 
z=3-4 for the FPL model. When we compare the dif- 
ference in total numbers of ejected BHs between Kl 
and lower recoil K2 realizations (upper panels in fig- 
ure 18), both the FPL and PGA model predict the 
the Kl case ejects up to 50% more black holes. The 
average masses of ejected black holes are similar in 
the FPL and PGA models (bottom panels in figure 
18) with the most common ejected black hole of 100 - 
1,000 Mq. 

The mass ratio of merging black holes sets the 
kick amplitude, with large differences between the in- 
coming black hole mass creating a small kick. The typ- 
ical number of ejected black holes per realization as a 
function of mass ratio is presented in figure 19. The 
FPL and PGA distributions are again very similar re- 
gardless of the kick distribution. The main difference 
is that higher mass ratio encounters are more com- 
mon in the PGA Kl model. Since the growth of black 
holes in the FPL model is very efficient - at Edding- 
ton rate - the incoming black holes grow faster and 
this decreases the mass ratio when the merger occurs. 
Hence, there are more mergers at q < 0.3 in FPL Kl 
and more mergers at q > 0.3 in PGA Kl. 

Some of the ejected massive black holes might 
contribute to the rogue black hole population. Micic 
et al. 2006 has shown that high redshift ejections of 
black holes from their host halos in the Local Group 
type environment could lead to their capture by Milky 
Way potential. Tracking these black holes are beyond 
the scope of this paper and will be addressed in future 
work. 



5 DISCUSSION AND CONCLUSIONS 

Using a very high resolution small volume cosmolog- 
ical N-body simulation (Via Lactea II), we construct 



150 



100 



50 



0, 



FPLKl 
PGAKl 
FPLK2 
PGA K2 



J L_ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 19. Number of ejected black holes per realization 
as a function of the mass ratio of the merging black holes. 
Black lines represent the FPL model, and red lines rep- 
resent the PGA model. Thick likes show the Kl distribu- 
tion, while thin dashed lines present the lower recoil K2. 
Since the Kl model yields on average an order of magni- 
tude larger kick velocities than K2 lower recoil model, the 
number of ejected black holes is larger for Kl. 



the merger history of black holes in low mass dark 
matter halos in a Local Group analog. We found that 
the method used to estimate the dynamical friction 
timescale makes a large difference in the total rate 
and redshift distribution of black hole mergers (see 
also HoUey-Bockelmann et al. 2010). In particular, for 
example, an N-body based dynamical friction estimate 
(Boylan-Kolchin 2008) yields longer merger timescales 
than does Chandrasekhar dynamical friction. These 
larger, more realistic merger timescales will postpone 
black hole mergers to lower redshifts. Moreover, if 
most of the black hole growth is tied to gas accre- 
tion that is activated by major mergers, then using 
Chandrasekhar dynamical friction will overestimate 
the black hole mass. 

We studied the growth and merger rate of black 
holes from seeds at redshift ~ 20 to low mass SMBHs 
at redshift zero, analytically incorporating the impor- 
tant sub-resolution physics to model the black hole 
dynamics and gas accretion. We find that having a 
constant mass black hole seeds, or a narrow BHMF 
leads to nearly equal mass ratio black hole merg- 
ers. These mergers have largest gravitational wave re- 
coil velocities, which eject the black holes from their 
host galaxies and suppresses SMBH growth. A wider 
range BHMF presents far fewer problems in growing 
SMBHs, even when the black holes spins are large. 

Many groups are beginning to explore the forma- 
tion and evolution of massive supermassive black holes 
in large volume cosmological volumes, with a sophis- 
ticated and self-consistent treatment of the sub-grid 
physics (e.g. Croton et al. 2006, Sijacki et al. 2007, di 
Matteo et al. 2008, Sijacki et al. 2009). With larger 
simulation volumes, these simulations are capable of 
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forming hundreds of massive SMBHs at the centers 
of massive eUiptical galaxies and clusters of galaxies. 
In fact, the FPL model we adopt is a semi-analytic 
model (Somerville et al. 2008) resulting from these 
techniques. Due to high computational cost, however, 
the same numerical approach has not been applied in 
smaller cosmological volumes where the bulk of the 
halo evolution is non-linear. Here, high redshift dwarf 
galaxies merge to form the Local Group galaxies, and 
massive black holes have masses well below lO'^ Mq . 

We use small volume nbody simulations and ap- 
ply two semi-analytic recipes for SMBH growth. One 
class of models is always limited by the black hole 
fundamental plane (FPL). Here, the growing black 
hole mass is a function of the total stellar mass in 
the galaxy bulge. This approach has successfully pro- 
duced SMBHs in massive galaxies. We show that the 
FPL model can also be used for black hole growth 
in small volumes. In particular, black hole mergers 
combined with a) gas accretion in a quasar mode and 
b) a time averaged Bondi-Hoylc accretion rate of 7 
xlO~* Meyr^^ in a radio mode, can produce the ob- 
served mass of the Sgr A* SMBH. Today, the observed 
Sgr A* accretion rate is an order of magnitude less ef- 
ficient than this radio mode accretion rate, but this 
may not be unreasonable, since the Milky Way is not 
considered an AGN anymore. Other low mass AGNs 
in the local Universe are observed to have 10 - 100 
times more efficient gas accretion. The FPL model, 
by design, guarantees the formation of a fundamental 
plane black hole. However, in reality, not every galaxy 
has such a system, due to the competition between 
black hole growth and nuclear star cluster formation 
in low mass spheroids. 

We have proposed an alternative method which 
could provide deeper insight into the formation of the 
lightest supcrmassivc black holes (<10''Mq) by incor- 
porating growth prescriptions calculated in small scale 
merger simulations. Gas accretion onto the black hole 
is always suppressed by stellar or AGN feedback, and 
depends on the structure of the gas, as well as the mi- 
crophysics (such as turbulence, fully rclativistic mag- 
netohydrodynamics, etc) that is currently physically 
beyond modeling. Our model provides a framework 
into which the microphysics can easily be incorporated 
once simulations and observations reveal the nature of 
these processes. In contrast, in the FPL model, these 
effects are irrelevant, since gas accretion occurs in the 
same manner in all < 10:1 galaxy mergers. The fi- 
nal outcome is always a fundamental plane black hole 
since feedback rapidly drives the evolution of funda- 
mental plane rather than the growth of black holes. In 
the PGA model, we can incorporate an accretion ef- 
ficiency parameter into the gas accretion prescription 
to make black hole growth less efficient for stronger 
feedback and vice versa. We use a formalism that re- 
lates the gas accretion efficiency to galaxy merger ra- 
tios; gas accretion is most efficient for major mergers 
and negligible for minor mergers (Cox et al. 2008). 



In addition, minor mergers will last much longer than 
equal mass mergers. When these two effects combine, 
the PGA model produces short bursts of high rate 
gas accretion for equal mass mergers, and long, slow 
gas accretion for 10:1 mass ratio mergers. As a result, 
most black holes in the PGA model never reach the 
Eddington rate. 

We also modeled the effect of gravitational wave 
recoil on assembling the lightest SMBHs. One way to 
mitigate gravitational wave recoil is to increase the 
black hole merger symmetry; smooth dense gas. for 
example, may quickly align the spin vectors of each 
SMBH before coalescence (Bogdanovic et al. 2007, 
Dotti et al. 2010). In order to understand how damag- 
ing gravitational wave recoil can be, and to pin down 
the epoch that a growing black hole is most vulnerable 
to recoil, we do not pre-select the merger configura- 
tions with the lowest recoils in Kl and select spins 
aligned to the orbital angular momentum in lower re- 
coil model K2. In previous work, when black holes are 
merging in massive stellar systems, the large gravita- 
tional potential creates large escape velocities. Hence, 
a K2 kick distribution with much lower recoil veloci- 
ties helps retain the central black holes (Volonteri et 
al. 2010). In our merger tree, most of the dark matter 
halos have masses below 10^" Mq, and the typical halo 
mass is 10^ Mq - 10^ M©. Since the escape velocities 
in these systems are below 100 kms~^, even the mod- 
erate kicks will eject the SMBH. When gravitational 
recoil is included in the FPL model, the final black 
hole mass is unaffected in most cases. This is the con- 
sequence of the FPL requirement of a very high gas 
accretion efficiency. Even in the most damaging Kl 
realization where the black hole is ejected twice from 
the merger tree, the accretion of gas is so efficient that 
a ~ 100 Mq black hole that is ushered to the center 
has time to grow to 2x10^ M©. When gravitational 
recoil is included in the PGA model, the final mass of 
the central black hole in the Kl realizations is lower 
by roughly an order of magnitude than in K2 (lower 
recoil model). 90 % of the K2 realizations are in the 
Sgr A* mass range. PGA model predicts that kicks 
have to follow lower distribution of gravitational re- 
coil velocities, K2. 

When we look at the mass of the black holes at 
the centers of local field dwarf galaxies, none of the 
models, FPL or PGA, matches the extrapolation of 
the Ferrarese (2002), relation which relates black hole 
masses to the masses of host dark matter halos. Even 
when we remove gas accretion as a mechanism for 
black hole growth, their masses will be larger than 
predicted by Ferrarese relation as the result of seed 
black hole mergers, despite the fact that the initial 
black hole seeds were Pop HI remnants and therefore 
considered light seeds. It is possible that the black hole 
fundamental plane exists for low mass SMBHs in field 
dwarf galaxies, but it may simply have a difi^erent slope 
than for high mass SMBHs; indeed this was alluded to 
in the original Ferrarese work. We do find a shallower 
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slope in all models. In FPL, the slope is close linear 
and tracks the following proportion: Mdmh / Mbh 
= 10^ The PGA model predicts an even shallower 
slope - 0.64 - 0.80 - although the scatter caused by 
gravitational wave recoil makes it hard to distinguish 
between FPL and PA. Note that gas consumption in 
field dwarf galaxies may favor nuclear star cluster for- 
mation rather than a massive black hole. Both FPL 
and PGA models are consistent with this picture since 
roughly half the dwarfs contain no central black hole, 
though the PGA models predict more empty halos on 
average than FPL. It is important to note that the 
black hole demographics in the dwarf galaxy regime 
is critical to distinguish between the FPL and PGA 
models. 

We also identified two massive rogue black hole 
populations. First, we find ejected black holes, or black 
holes with kick velocities larger than the host halo es- 
cape velocity. We will address the properties of this 
population with a detailed treatment of their dynam- 
ics in the future. The second arc rogue black holes from 
the remnants of satellites that had merger timcscalcs 
longer than a Hubble time. We find hundreds of these 
rogue black holes at distances of ~ 200 kpc from the 
galactic center and with masses from 10- 10^ M© (de- 
pending on the model). The most massive of these may 
be detected in X-ray and are possible mesolensing can- 
didates. These two populations of black holes, which 
have quite different formation mechanisms, might not 
be easily disentangled. It is possible that a growing 
primary halo may recapture black hole ejected by re- 
coil, in which case both populations occupy a similar 
parameter space (Micic et al. 2006). 

In conclusion, our model successfully reproduces 
Sgr A* without limiting the black hole growth to the 
black hole fundamental plane. It predicts that kicks 
have to follow lower recoil velocity distribution (K2). 
It also predicts a substantially fiatter mass function 
(slope = 0.64 0.8) for lower mass SMBH in local field 
dwarfs. It predicts Milky Way rogue black holes, and 
a population of ejected black holes. This model pro- 
vides cosmological framework into which the results 
of small scale simulations can easily be incorporated 
in the future. At this point, the model is limited only 
by the lack of understanding of how various feedback 
mechanisms couple with each other and with the gas 
accretion onto the black hole. There is ample ground 
to cover in this respect. As one example, Ostriker et 
al. (2010) claim that all previous studic^s of SMBH 
growth may overestimate the final black hole mass by 
~ two orders of magnitude by treating only the energy 
part of AGN feedback, neglecting the importance of 
momentum and mechanical feedback. Hopefully, this 
will be resolved in the future. 

One issue is how well the VL-2 simulation rep- 
resents the Milky Way merger history. The reality is 
that VL-2 represents one out many possible evolution- 
ary histories of the Local Group, and not necessarily 
the correct one. The fact that only major mergers acti- 



vate the quasar mode means that the number of major 
mergers each galaxy goes through has a huge impact 
on the black hole accretion history. This underlines 
the importance of approaching the problem of Sgr A* 
growth statistically by simulating hundreds of highly 
resolved small simulation volumes. To resolve thou- 
sands of Milky Way mass halos with VL-2 resolution, 
gas physics, and black hole dynamics implemented on 
the fly will be computationally too expensive in the 
near future. 
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